Casimir force due to condensed vortices in a plane 
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The Casimir force between parallel lines in a theory describing condensed vortices in a plane is 
determined. We make use of the relation between a Chern-Simons-Higgs model and its dualized 
version, which is expressed in terms of a dual gauge field and a vortex field. The dual model can 
have a phase of condensed vortices and, in this phase, there is a mapping to a model of two non- 
interacting massive scalar fields from which the Casimir force can readily be obtained. The choice 
of boundary conditions required for the mapped scalar fields and their association with those for 
the vectorial field and the issues involved are discussed. We also briefly discuss the implications of 
-^Vj- our results for experiments related to the Casimir effect when vortices can be present. 

PACS numbers: 03.70.+k, 11.15.Yc, ll.10.Ef 

C/3 1 I. INTRODUCTION 

The Casimir effect is a manifestation of the quantum vacuum fluctuations that can be tested at mesoscopic scales. 
This quantum phenomenon has been of interest to fundamental physics since its prediction by Casimir in 1948 1] 
and it has been studied extensively both theoretically and experimentally since then (for a recent review, see e.g. 
Ref. 0). In particular, many experiments have been measuring the Casimir force with increasing precision. It is 
^ • then of increasing interest to look for possible novel situations where the theoretically computed Casimir force can 
be confronted to experiments and in which related quantum phenomena, associated with the quantum vacuum, can 
w—. then be probed and tested in the laboratory. 

In this work we want to study how a vacuum state made of topological excitations, more precisely, a vacuum 
. constituted of condensed vortices, will affect the Casimir force between perfectly conducting parallel lines in a plane. 
1 Let us recall that stable vortex configurations can appear in important condensed matter systems, like in high- 
\& . temperature superconductors and superfluids (for a detailed presentation see e.g. Ref. Q and references therein). It 
is expected that the nontrivial vacuum state, where vortex condensation may be present, to be of relevance in many 
recent experiments as far the Casimir effect is concerned. For example, it has been pointed out in Ref. [3| that unusual 
. ■ behaviors of superconductors may be found when the size of the samples shrink. But we note that this is precisely 
0^ the case, at the nanometer scales, that we expect that the Casimir force to become more appreciable. In the case 
of vortex-based superconducting detectors 0-0] for instance, it can be expected that the Casimir effect can possibly 
alter the microscopic parameters of the detector, analogously to the case reported in Ref. 

The simplest and, in our opinion, the most direct way for studying a vortex condensate state is through the use 
of dual transformations involving the field variables of the original Lagrangian density. Dual transformation is a 
common technique through which the topological excitation content, which in our case we will be concerned with 
vortex excitations, of a given theory can be made explicit. This is reminiscent of similar approaches first used in 
condensed matter studies performed on the lattice [9J and of routine use since then. Through a series of appropriate 
dual transformations involving the original fields in the functional action, an equivalent action is obtained where the 
vortex excitations are made explicit. By properly matching our dual action to a field theory model, it is then possible 
to write it in terms of a vortex field coupled to a gauge field (for earlier implementations of this procedure, see for 
example the work done in the Refs. [R>ll2j and references therein). 

Here we will show that the calculation of the Casimir force in a vortex condensed vacuum will be equivalent 
to computing the force in the case of a massive gauge field in a Maxwell-Proca-Chern-Simons (MPCS) analogous 
model. As shown in [HI, [Hj], the MPCS model can be seen as the dualized version of a Chern-Simons-Higgs (CSH) 
model, where the vortex excitations of the CSH model are made explicit. Vortex condensation in Chern-Simons (CS) 
type theories, particularly in self-dual models, has been shown possible for some critical value of the Chern-Simons 
parameter (l3.ll4|. with the determination of the condensation point explicitly obtained in 12]. The Casimir force for 
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the dual MPCS type of model is then here studied in this context, which is valid deep inside the vortex condensate 
phase. 

Irrespective of the connection between the dual MPCS model with the CSH model, this study has an interest of 
its own associated with the determination of the Casimir force for massive vectorial fields. Recall that the MPCS 
model represents, by itself, massive photons in 2+1 dimensions, with the photon mass having contributions from 
the usual Proca and the CS terms. While the mass coming from the Proca term can be seen as been generated 
though a symmetry broken scalar field, the one from the Chern-Simons term is of purely topological origin [l5| . 
The issue of the Casimir force for a vectorial field is closely related to important questions, both from experimental 
and theoretical points of views. For instance, in the case of 3+1 dimensions, the authors of Ref. |16[ have analyzed 
the existence of new expressions for the electromagnetic field between conducting plates, where the photon has a 
possibly non null mass, and study the dependence of the Casimir force on this mass. Later, in Ref. [17] . the Proca 
equations where used to represent the photon mass. In Ref. Q, it was considered the mass acquired by the photon 
due to the spontaneous symmetry breaking that takes place when a superconducting detector passes from its normal 
(N) to the superconducting (S) state, as a consequence of the detection of an external photon. In that reference, 
it was argued that the Casimir effect can alter the S-N transition in a detectable way and altering in this manner 
microscopic parameters of the detector. It has been stated by the authors in Q that the standard expressions for 
the electromagnetic Casimir force (considering massless photons), the Maxwell equations should be replaced by the 
Proca relations for massive photons. 

The Casimir force between perfectly conducting parallel lines in a plane for a MCS theory has been determined 
previously in [T§| - [2l| . In particular, it has already been shown in Ref. [lH that the Casimir force obtained is identical 
to that coming from a massive non-interacting scalar field in 2+1-dimensions 22]. This result can be understood from 
the fact that in both theories the field satisfies the Klein-Gordon equation of motion and both have only one degree 
of freedom. Note that in this case, the Chern-Simons term provides a mass term for the gauge field, but this is a 
topological mass that still maintains the field with only one (transverse) polarization degree of freedom. Besides, the 
boundary conditions (BC) in both models can be matched, making the Casimir force in both models to agree. On the 
other hand, in a symmetry broken case, a Proca mass term is generated for the gauge field and this in turn acquires 
a longitudinal polarization degree of freedom, in addition to the transverse polarization. For the case of the MPCS 
theory, the gauge field now has two polarization degrees of freedom. Similarly to the case of the MCS theory, we now 
expect that the respective Casimir force could be related to two massive noninteracting scalar fields. In fact, it is well 
known in this case, that the quantum mechanical analogue of the MPCS theory is equivalent to two noninteracting 
harmonic oscillators with distinct frequencies [ll|. At the quantum field theory level this must then correspond to 
the case of two noninteracting massive scalar fields. This has been shown explicitly in [23| . where a mapping between 
the two theories was constructed. However, in order to associate the corresponding Casimir forces for both theories, 
a careful consideration of the boundary conditions must be accounted for. This issue was earlier discussed in Refs. 
[Ill, [24]]. Here we will give a detailed account for the issue of the boundary conditions when mapping our dual MPCS 
theory with vortices with that of the two noninteracting massive scalar fields mapped model. This will allow us then 
to readily obtain the Casimir force for the model we have studied here. 

The remaining of this work is organized as follows. In Sec. [Hi we introduce the MPCS model as a particular limit 
of the dual vortex model considered in Ref. [l2j and summarize the relevant equations and relations that will be of 
relevance for this work. In Sec. IIII1 we discuss the mapping that leads from the initial MPCS theory to a model of two 
massive and noninteracting scalar fields. We also analyze the respective mapping between the boundary conditions 
needed for those two models. In Sec. IIVI we then derive the Casimir force related to a vacuum state of condensed 
vortex excitations from the dual MPCS theory considered and contrast the result with the case of absence of vortex 
excitations. In Sec. [V] we give our concluding remarks and discuss possible extensions of our work. One Appendix is 
included to show some of the technical details. 



II. THE MPCS THEORY AS A DUAL MODEL FOR VORTICES IN THE PLANE 



Let us initially consider the CSH model in 2+1 dimensions, written in terms of a complex scalar field and an Abclian 
gauge field, which here we will represent here by 77 and h^, respectively. The quantum partition function and the 
action of the model have the forms (in Euclidean space-time and with indices running from 1 to 3) 



Z = J Dh^Dri'Dr]* exp {-S E [h^ 77, r)*}} 



d 3 x 



-i-e^h^ + \D^\ 2 + V{\r,\) 



(2.1) 
(2.2) 
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with = dfj,h u — d v h^, = + ieh^ and 6 is the CS parameter. V(M) i s a symmetry breaking polynomial 
potential for the scalar field, independent of its phase, which has a non-null vacuum expectation value (VEV) \(t))\ = 

v 7^ 0. As examples of V^M)' we can c ^ e tne usual quartic potential V(|7y|) = A (|ry| 2 — v 2 ) 2 /4 and the sixth-order 

self-dual potential [IBj]: V(M) = e 4 (|?7| 2 — v 2 ) 2 1 77 1 2 / O 2 . Writing 77 in a polar form rj = (p/V2) exp (ix), the VEV for 

r) becomes v = po / 

The field equations for and rj are known to have a nontrivial solution associated with a vortex field configu- 
ration [2^|. When expressed in polar coordinates (r,x), the nontrivial solution can be put in the generic form that 
represents charged vortices: 



Wtcx = £,{r) exp(mx) , (2.3) 

Ti 

Vvortcx = -H r ) 9aX ) (2-4) 

e 

where n is an integer that can be interpreted as the vortex topological charge (or winding number) and (£(r), h(r)) 
are obtained by numerically solving the classical field differential equations, subjected to the BC: 



lim £(r) = 0, lim £(r) = v, (2.5) 
lim h(r) = 0, lim h(r) = 1. (2.6) 

i — s-0 r—too 

A vortex represented by Eqs. (I2.3|) and (|2.4I) can be seen as carrying an "electric" charge Q (the spatial integral of 
the 0-component of the density current j M ) attached to a magnetic flux $ given by 



$ = j d 2 xH 12 = | . (2.7) 

This fact is a direct consequence of the presence of the Chern-Simons term and also implies in an anyonic behavior 
of the charge-flux composite, which has spin s = Q$/(47r) 26]. It can be also demonstrated that when r approaches 
to infinity (or sufficiently far from the vortex core), the flux $ becomes quantized, being given by an integer multiple 
of flux quantum (26[: $ = 2nn/e. 

The vortex degrees of freedom present in the original theory Eq. (|2.2|) can be made explicit through a series of 
duality transformations jTl], [l3 - The final result is a theory of the form of a Maxwell- Chern-Simons-Higgs (MCSH) 
model, where the vortex solutions, represented by Eqs. (|2 .3[) and (|2.4[) , are associated with particles represented by a 
complex scalar field ip that is coupled to a dual vector field A^. The original fields and the dual fields in the classical 
level are related to each other e.g. by p (d^x + e/i^) = (o / e)e tiV1 d ll A 1 / (2it), where a is an arbitrary parameter with 
mass dimension. The resulting dual Euclidean action becomes equivalent to a MCSH theory of the form [ll|, G3 : 



S, 



dual 



d 6 x 



16tt 



2 o2p 2 ^ 



F'L + i 



87r 2 e 



^awy -Audi/ A*y -\- 



2<7 

e 



(2.8) 



where — d^Ay — d v A^, V{\ip\) is the effective potential term for the vortex field and Lq is a gauge fixing term. 
Note that in the dual model, Eq. (I2.8[) . the new CS coefficient appears inversely proportional to the initial one in 
Eq. (|2.2p . G — > — l/(47r 2 6). This dualization of the CS coefficient is a consequence of the transformations used (see 
also Refs. [23, H). 

As argued in Refs. [HOJ], there is a critical value for the CS coefficient in the CSH theory below which vortex are 
expected to be energetically favorable to condense. In terms of the dual action (|2.8|) . this can be expressed in terms 
of an existence condition for a VEV for the dual vortex field, given in terms of the first derivative of the potential 
with respect to the vortex field, V v ' ortex (|V'| = ipo/v*) = 0, or, analogously, that the quadratic mass term in the vortex 
potential be negative below some critical C , with C determined by the condition on the second derivative of the 
effective vortex potential with respect to the vortex field, V v " rtex (0 = O c ) = 0. In Ref. [12J this critical value has been 
obtained as given by O c ~ (e 2 /7r)ln6 and shown to be robust against quantum corrections, that changes 9 C by no 
more than about 17%. In this work we are interested in deriving the Casimir force starting from the dual action (|2.8[) 
in the case where vortex condensation is favorable, i.e, for the region of parameters where < C . 

Since the Casimir force is associated with quantum vacuum fluctuation, if we want to determine an expression for 
that force in the case of the MCSH theory of the form of Eq. (|2.8|) . we can, as an approximation, consider only 
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small variations of the vortex field around its nontrivial constant VEV ipo- I n other words, if we are deep inside the 
vortex condensed phase, fluctuations of the vortex field can be neglected, much like in the London approximation in 
condensed matter problems Q. This approximation can then be seen as a limiting case for Eq. (|2.8[) . in which the 
term |<9 M V + 2i<rA^ip/e\ , in Eq. (|2.8p . gives rise to a Proca like term (i.e., we are in the vortex symmetry broken 
phase in the dual action), which can then be written simply as rr? A p A^. We can also make use of the arbitrariness 
of a in Eq. (|2.8p to rewrite the action in the form of a MPCS model (see also the Appendix). Considering a = 2irepo 
and going back to Minkowski space-time, the corresponding MPCS Lagrangian density then becomes, 



£ = -\f^F^ + l -m 2 A»A^ + ^ X A^A X . (2.9) 



where in Eq. (|2.9p . for convenience, we have redefined the parameters as 



H = 2e 2 po/©, m = 4 7 r / o o <0o . (2.10) 

The association of a covariant derivative of a field with a mass term for a boson, in the broken vacuum state (for 
the dual theory in our case), is a well known result, very similar to the mechanism of mass generation for photons 
inside superconductors, which can be explained in terms of a symmetry breaking like in the Landau-Ginzburg model 
for superconductivity [29(. 

In the next section, the Casimir force associated with the theory described by Eq. (|2.9|) is determined by noticing 
that it can be mapped to an equivalent model of two non-interacting massive scalar fields, as mentioned in the 
introduction, and by properly choosing the appropriate boundary conditions for the corresponding scalar fields and 
the dual vector field. In the association between the two theories, the two initial massive degrees of freedom of the 
MPCS model are transposed to two degrees of freedom represented by the two scalar fields, as it should be expected 
in the quantum mechanical analogue case [l5j . 



III. THE EQUIVALENT MODEL AND THE MAPPING BETWEEN THE BOUNDARY CONDITIONS 

The MPCS theory given by Eq. (|2.9[) can be mapped, after a sequence of mathematical transformations, in a model 
of two noninteracting massive scalar real fields (<f> and ip) in 2+1 dimensions |23l j. Next we will explain the main steps 
needed for the mapping and that will be useful for setting the respective BC needed in the calculation of the Casimir 
force. 



A. The MPCS theory equivalence to two noninteracting scalar fields 



From Eq. (|2.9[) . the Euler-Lagrange equation for the dual gauge field A^ is 



dpF" + m 2 A a + ^t apf) F pfi = 0, (3.1) 



while the canonical momenta are 



n o = d^ = Q » = J^L = F i0 + —e ij A v (3.2) 
dA° dA* 2 J ' y J 

where the indexes i and j vary from 1 to 2. The relation 7r° = is a primary constraint of the model, which also 
shows a secondary one, given by 



drf + Je ij >ij + m 2 A « 0. (3.3) 

The primary and the secondary constraints are directly related to the reduction of the number of degrees of freedom 
of the system (from 3 to 2). We also note that the secondary constraint permits that one can write Aq in terms of 
the components A+. This possibility can be seen as a direct consequence of the fact that the vectorial field mass m 
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is non-null. As a consequence of the constraints, the physical degrees of freedom of the system are represented by Ai 
and 7Tj. The quantum partition function can now be written in the form 



Z = J T>A l D-n 1 exp 
where is the physical Hamiltonian density, 



(3.4) 



where 2£y, Qij and S lJ are defined by 



(3.5) 



_ didj 

A ■■> ~ + ,,, J : 



(3.6) 



Qij — n \ ~\~ o didj 

2 \ m 2 



di = eijd 3 , 



(3.7) 



S ij 



Am 2 



[(v 2 - 



m 2 ) g ij + d l d j ] 



(3.8) 



It is important to note that in order to write the Hamiltonian density in the form of Eq. (|3.5[) . the surface terms 
generated by the integrals of di^djir^), di{ir l e> k Fjjf) , d l (A^diAj) and d l (A : >djAi) are neglected. As we will show 
later, this can be shown to be indeed the case for the boundary conditions that we will be considering here. 

Next, we introduce two new variables, A\ and fr* (i = 1, 2), defined by the relations [23| 



A x = 
A 2 = 

7T 1 = 



(Sr x ii - 6j l A 2 ) /(20), 
- (o^Ax + b^A 2 ) /2. 



where 



4 ' 



and Oi and 62 are operators whose squares are given by 



(3.9) 

(3.10) 
(3.11) 
(3.12) 



(3.13) 



of = ( --e*K xl 



>12 



s 2 



(3.14) 



6- 



1 



-e 2 K n + 6Q 12 + S 



2'2 



(3.15) 



We note from the above equations that when acting the operators 61 and 2 on some function (e.g. 4>(x)), they 
cannot simply be written in terms of the derivatives of the function. In Eqs. (|3.9[) and (|3.12l) . fr* and Ai can be seen 
as intermediate variables, related to the fields {0, <p} and their respective momenta {n^, as 
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7Ti 


1 


v 2 ; 




= 71^" 




M 


= V24,, 




A 2 


= V2<p. 





(3.16) 

V, (3-17) 

(3.18) 
(3.19) 

The set of mathematical transformations shown above make it possible to rewrite the Hamiltonian of the MPCS 
model as a sum of two separated and independent Hamiltonians associated with two noninteracting scalar fields tp 
and 4> like 

K = \ [4 + <K™? - V 2 )0] +IK + v(m 2 2 - V 2 )^] , (3.20) 

where 

m 1 = 0-£ m 2 =d+^. (3.21) 

The relation between the model described by Di, Eq. (I3.20[) and the MPCS theory, can now be used to obtain the 
Casimir force for the dual model Eq. (|2.9j) . the analogous model describing condensed vortices in the dual formalism. 
Since the Casimir force for a massive scalar field in 2+1-dimension is well known [22|, we can use this mapping between 
Eqs. (|2.9I) and (|3.20p to readily obtain the Casimir force from (|2.9I) . provided well defined BC are considered. Thus, 
we next turn our attention to this issue of the BC for the mapped theory. 

B. The BC mapping between the gauge held and the scalar fields 

The problem of the BC is a non trivial problem for the case we are considering here. This is because some usual 
BC considered for scalar fields cannot be directly written in terms of the vectorial field j4„, considering the form of 
the used mathematical transformations. This is an important issue, since it is well known that the Casimir force (for 
both its modulus and orientation) depends significantly on the BC considered. 

To set the BC that will be considered for the vectorial field, we must map mathematically the set {<f>, ip} (and/or 
the derivatives of those fields) into {Ao, A\,A 2 } (and/or the derivatives of those fields components). Hence, we need 
to invert the relations {Aq, A±, A 2 } — > {4>, ip} given in the previous subsection. First, it is easy to see that 



A l = 



61 V ~ 62 Vl /(V20). (3.22) 



In the momentum space, the expression (|3.22p can be easily understood but, since the physical BC are not taken in 
the momentum space and our aim is to relate the BC considered for the scalars and the vectorial field components, 
we must first of all understand the meaning of OfV and O^ 1 ip in the real space-time (at least at the boundaries) . In 
order to do this, we initially write the explicit forms for the operators K\\, Q\ 2 and S 22 that appear in the expression 
for O 2 , Eq. (|3.14[) . After that, we invert both sides of that equation in order to find an explicit expression (in terms 
of derivatives) for Oj~ 2 0. We can view this operation as a step to find a clear interpretation of the action of O^ 1 in 
(j), a necessary task if we want to write Ai in terms of 4> and V 5 using Eq. (|3.22l) . Thus, we can write 

6^ =(A- Bdf) 0, (3.23) 
where we have used the definitions of two new constants: 

T~T + ~l + 4^i' ~2^~2^^2 _ 4^- {3 ' 
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If we write 1 2 as the square of a certain operator, then it is possible to identify that operator with . By this 
manner, we can have an expression for O^ 1 ^, at least at the boundaries. With this aim, we will consider that the 
system is constrained in a infinite strip, with boundaries at x = and x = a. Also, we will consider that the field <f> 
obey the Neumann BC, d\<j)(x = 0) = d\<f>{x = a) = 0. So, at the boundaries, we can rewrite Eq. (|3.23p as 




(3.25) 



and hence 




Thus, it is now easy to write the first term of the LHS of Eq. Q3.22p at the boundaries, in the real space. For the 
second term in the LHS of that equation, we use the same arguments. Considering Neumann BC for <p, we can write, 
at the boundaries, 



6 2 - V = Vctp, 



(3.27) 



with 



6 2 Oil m 2 , 



A) 

4m 2 I 



(3.28) 

/(V26) and so we see that A± 



From our results, we can write Eq. (|3.22|) . at the boundaries, as A\ — 

must also obey Neumann BC. 

Making analogous considerations for A 2 and making use of the Eqs. (|3.10p . (13. 16)) and (|3.17[) . we obtain that 



A 2 = 



(6l7T + 6 2 7T V 

71 



\/20i 



Q22 
2 



■-V2O2 



Q22 
2 



tp. 



We can now use the Eq. (|3.29l) to determine the behavior of A 2 at the boundaries, where 0i</> 
Noticing that we are considering Neumann BC for <p and ip, using the Hamilton equations {71^ 
writing the explicit forms of the operators Q22 and S* 12 , we conclude that at the boundaries 



(3.29) 



'A 



, 2 if = tp/VC. 
Kip = dotp} and 



A 2 = 



d 4> d <p 



>2A 



l2C 



(3.30) 



The Eq. (|3.30[) implies that A 2 must also obey Neumann BC. 

The BC considered for A\ and A 2 , together with the Euler-Lagrange equations and the definitions of the canonical 
momenta, define the components of the strength tensor at the boundaries. The behavior of those components of 
should not be confused with new BC imposed on the vectorial field, but just another way to write the BC themselves. 
For instance, from the definition of 7r\ we have 



F 20 = vr 2 + 



(3.31) 



or yet, from Eqs. f3~3]h (|5T2"j) . (j3~T5|) and 



F 



20 = V A _ 



6rV 



(3.32) 



So, at the boundaries, we have 
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F 



V2 
2 



Acj) + VCip 



(3.33) 



which implies that diF 20 = at the boundaries. 

It is worth to note that we can also write that condition in terms of the dual tensor F^, 



2 vp ' 



(3.34) 



in the form 



d 1 F 1 (x = 0) = dxF 1 (x = a) = 0. (3.35) 

Let us now investigate the implications that the Neumann BC for Ai, A2 and F 1 have on the behaviors of Ao, F° 
and F 2 at the boundaries. First, writing Eq. Q3.1j) for a = and a = 1 and using the BC already found, we get 

(d^ 1 + m 2 ) A - ^d 2 A 1 + d 2 F 20 = 0, (3.36) 
(d d° + d 2 d 2 + m 2 ) A 1 - daSp-Ao + iiF x = 0. (3.37) 

From Eq. Q3.37p . it is easy to see that dodid 1 A = 0. If we represent A in terms of its Fourier transform in y and t 

M*, V,t) = J f^ - *"* J 7^ eihu Mx, k > w), (3-38) 

then we easily obtain that dodid 1 Ao = (valid for any t), which implies that did 1 Ao = at x = and x — a. We 
can now use this relation in Eq. (|3.36l) to obtain 

m 2 A - fJ ,d 2 A 1 + d 2 F 20 = 0. (3.39) 

Since A\ and F 20 are subjected to Neumann BC, Eq. Q3.39[) implies that Ao is also subjected to the same BC as well. 

By considering that F° is subjected to the requirement doF° = [l8j], then the Bianchi identity d^F^ = and 
the Neumann BC imposed on F 1 imply that dod 2 Ai = at the boundaries. We then conclude, after straightforward 
considerations, that A\ must be constant at the boundaries. Using Eqs. (|3. 221) . (|3.26p and (I3.27[) , we can write, for 
x = and x = a, that 




(3.40) 



By noticing that <p and ip are independent functions and that A\ is constant at the boundaries, we can conclude, using 
Eq. p. 401) . that 4> an d f must be constants at the boundaries as well (their values cannot depend on the explicit 
values of x, y, t). This fact has direct implications for the behavior of F M , as we next show. First, using Eqs. (|3.40[) 
and f|3 . 33[) . we can write 



F 20 {x = 0) = F 2Q (x = a) = -1- (yfAmit + VCm 2 <p) . (3.41) 

We see then that F 1 must be constant at the boundaries. Using Eq. (|3.34l) and remembering that A 2 and Aq obey 
Neumann BC, we conclude that F° and F 2 are null at the boundaries. An important particular case of the BC to 
be obeyed by the components of the dual tensor is F° = F 1 = F 2 = 0, which characterizes perfect conductor BC. 
In this case, using the definitions of the canonical momenta and the considerations about the behavior of Ai at the 
boundaries, it is easy to prove that the surface terms generated by the integrals of d^^djir^), di(n t e' k Fjk) ) d l (AjdiAj) 
and d % (A J djAi), that appears in the generating functional, will give no contributions. This justifies neglecting these 
contributions, as we have assumed earlier. Analogously, the perfect conductor BC used here allow us to neglect the 
surface terms involving (f>, (p and their respective conjugate momenta, so to obtain the final Hamiltonian density 
Eq. 
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IV. THE CASIMIR FORCE 



Once we have fixed the required BC, we can now easily find the Casimir force for the dual MPCS theory Eq. (|2.9|) . 
This follows directly from the equivalence between the original theory Eq. (|2.9|) with the model represented by 
Eq. (|3.20p . The Casimir force for a massive scalar field subjected to Neumann (or Dirichlet) BC in 2+1-dimensions 
(which is also the same as the one computed for a MCS theory) is [H, [H[ 

1 f°° y 2 

/scalar (m s , a) = — 3 / dy — , (4.1) 

16ira 6 J 2msa e y -1 

where m s is the mass of the scalar field. Using the equivalence between Eq. (|2.9|) with Eq. p. 201) . we can then 
immediately write the corresponding Casimir force, in the presence of a vortex condensate, as 

/vortex = /scalar (^1, a) + /scalar (™2 , d) , (4.2) 

where mi and m2, using Eqs. (|2.10p . (13. 13)) and (|3.2ip . are given by 



e 2 pl ( F 167T 2 ^6 2 \ 
For small values of mass, ma < 1, Eq. (|4.1I) can be expressed as 



/scalar (m s , a) 



87ra 3 



C(3) - (am s ) 2 + 



2(am s ) 3 (am s ) 4 



5„,5x 



G 



0(a°TO 



(4.4) 



where £(z) is the Riemann zeta-function. Using Eq. (|4.3|) and keeping for simplicity up to the quadratic term in the 
mass in Eq. (|4.4I) . we obtain for the Casimir force Eq. (14.2[) the result 



Ana 3 



C(3) 



e 



eV 



(4.5) 



The result (|4.2I) allows us to immediately conclude that in the presence of vortex matter (^ 7^ 0), the Casimir 
force is always smaller in magnitude than in the absence of vortices. In Fig. [T]we show the Casimir force Eq. 
as a function of the vortex condensate. 




This result can be interpreted as follows. Vortices are expected to repel each other, much like as in the standard 
mean- field phenomenology for type-II superconductors when vortices can form [30|, e.g. in the Shubnikov (vortex) 
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phase, where above some critical magnetic field vortices are present. The repelling vortices will exert an opposite, 
repulsive force on the external conducting lines that tend to counterbalance the attractive Casimir force, tending 
to make it smaller the larger is the VEV of the vortex condensate. The resulting Casimir force can then be made 
sufficiently small in the presence of vortex matter, though it will never be exactly zero or become repulsive, as it can 
be clear from the expression for the Casimir force and from Eq. (|4.3p . where of course m 1 ^ 2 ) > 0- 

V. CONCLUSIONS 

We have studied in this work how a nontrivial vacuum state, with condensed vortex excitations, affects the Casimir 
force between two conducting lines in a plane. By starting from a CSH theory with field equations having vortex 
solutions, and using its dualized form, which results to be analogous to a MCSH theory, vortex degrees of freedom 
are made explicit. In the vortex condensation regime, the dual model can be expressed simple as a MPCS like theory, 
which in turn can be mapped in a two noninteracting massive scalar field model. Using the known expression for the 
Casimir force for a massive scalar field, the corresponding Casimir force for the case of vortex matter between the two 
lines has been computed. 

We have shown that the Casimir force in the presence of vortex matter is smaller than in the absence of vortices. 
This result may have implications for Casimir force experiments using e.g. superconductors, like in the next generation 
of experiments [3lj . In particular, the presence of vortices can make the Casimir effect to be much smaller, making 
its detection through measurements possibly more difficult. 

Another important issue that could be considered is the possibility of using our results to find the Casimir force 
for the MPCS theory in the case of moving boundaries (i.e., the dynamical Casimir effect). As mentioned in the 
introduction, it is expected the Casimir effect to play an important role in superconductors, specially at nanometer 
scales. Recently, the first experimental observation of the dynamical Casimir effect in a superconductor circuit (32j 
has brought great attention to this matter. Some of the considerations that we have used here are also valid in the 
dynamical case. Of course, how we set the boundaries (e.g. at x — and x — a) arc of decisive importance for 
determining the expression for the Casimir force. However, the mapping between the initial MPCS theory and the 
model of the scalar fields makes use only of the values of the derivatives of the functions <f>, ip and Ai at the boundaries. 
But the value of x itself at those boundaries is never actually needed at any step. In other words, the mapping used 
here is expected also to be valid in the case of moving boundaries, as long as the BC remain valid (e.g. for a perfect 
conductor parallel to the y axis, with movement in the x direction). Hence, we conclude that we can, in principle, use 
the same arguments we have considered to also study the dynamical Casimir effect for the MPCS model. However, 
to find the Casimir force in that case, we must know the force for a massive scalar field between moving boundaries, 
which is an issue that we intend to treat in a future work. 



Appendix A: The Energy-Momentum tensor and the Casimir force 

The Casimir force for the dual theory is expressed, as usual, in terms of the VEV of the T 11 -component of the 
symmetrized energy-momentum tensor [18l |2lj: force/length = (0|T n |0). Thus, we can first write 

T" v = F^F" + m 2 A ti A v - 'Lf u (F X F X + m 2 A a A a ) , (Al) 

where, for the sake of simplicity, we made use of the definition of the dual tensor F M , F^ = ^e^ vp F vp . Usually, the 
components F M are associated with the components of the "electric" and "magnetic" fields (F 1 = —E y , F 2 = E x 
and the scalar B — F°). In this work, E x , E y and B may or may not (in the case of the dual gauge field) represent 
a physical massive electromagnetic field (we are just borrowing the usual nomenclature). Hence, ^0|T n |0) at the 
boundaries can be written in terms of VEVs of products like A^ A v and derivatives of them, taken at x — or x = a 
(the explicit values of x at the boundaries will not be necessary for our purposes). Following [3], we write those 
VEVs, at x = 0, as 

(0\Ai t (x)A lf (x)\0) | X1=0 = Um {0\A»(x)A»(x')\0) , (A2) 

xi— ^x^— 

where x and x' stand for points in the three-dimensional space-time. But the VEVs in the RHS of Eq. (| A2|) are 
the two-point functions of the model, which can be written in terms of the functional derivatives of the normalized 
generating functional Z[J a ], where J a is a source: 
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(0\A"(x)A v (x')\0) = - 



6*Z[J} 



SJ^x) 5J v {x') 



J a =o 



SJ^x) 5J v {x') 



J DA/s exp (iS + iJ J a A° 
J DA 7 exp{iS) 



(A3) 



J a =o 



The derivatives appearing in Eq. (IA3|) are independent of the norm of the field. This shows that the Casimir force 
should also be independent of the normalization of the fields. For example, the arbitrary mass parameter a appearing 
in Eq. (|2.8|) can be put as a global multiplicative constant a 1 in all terms in Eq. (|2.9[) and be re-absorbed in a 
redefinition of the norm of A u . 



Acknowledgments 

C.R.M.S. was supported by Coordenagao de Aperfeicoamento de Pessoal de Nfvel Superior (CAPES-Brazil). 
R. O. R. is partially supported by research grants from Conselho Nacional de Desenvolvimento Cientffico e Tec- 
nologico (CNPq) and Fundagao Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro (FAPERJ). 
J. F.M.N, thanks J. A. Helayel Neto and D. T. Alves for very fruitful discussions on the subject of this paper. 



[1] H. B. G. Casimir, Indag. Math. 10 (1948) 261 [Kon. Ned. Akad. Wetensch. Proc. 51, 793 (1948)] [Front. Phys. 65, 342 

(1987)] [Kon. Ned. Akad. Wetensch. Proc. 100N3-4, 61 (1997)]. 
[2] M. Bordag, G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Advances in the Casimir Effect (Oxford University 

Press, Oxford, 2009). 

[3] H. Kleinert, Gauge fields in Condensed Matter, vol. I: Superflow and Vortex lines, (World Scientific, Singapore, 1989); H. 

Kleinert, Gauge fields in Condensed Matter, vol. II: Stresses and Defects, Dif- ferential geometry, Crystal defects, (World 

Scientific, Singapore, 1989). 
[4] Alan T. Dorsey, Nature 408, 783 (2000). 

[5] A. M. Kadin, M. Leung, A. D. Smith, J. M. Murduck, Appl. Phys. Lett. 57, 2487 (1990). 
[6] A. D. Semenov, et. ai, Supercond. Sci. Technol. 20, 919 (2007). 
[7] A. D. Semenov, et. ai, Physica, C 468, 627 (2008). 

[8] D. Brandt, G. W. Fraser, D. J. Raine and C. Binns, J. of Low Temp. Phys. 151, 25 (2008). 
[9] M. E. Peskin, Annals Phys. 113, 122 (1978). R. Savit, Rev. Mod. Phys. 52, 453 (1980). 
[10] H. Kawai, Prog. Theor. Phys. 65, 351 (1981). S. -J. Rey, Phys. Rev. D 40, 3396 (1989). R. O. Ramos, J. F. Medeiros 
Neto, D. G. Barci and C. A. Linhares, Phys. Rev. D 72, 103524 (2005) [hep-th/0506052]. R. O. Ramos, D. G. Barci, 
C. A. Linhares and J. F. Medeiros Neto, Braz. J. Phys. 37, 251 (2007) [arXiv:0704.1117 [hep-th]]. 
[11] Y. Kim and K. -M. Lee, Phys. Rev. D 49, 2041 (1994) [hep-th/9211035]. 

[12] R. O. Ramos and J. F. Medeiros Neto, Phys. Lett. B 666, 496 (2008) [arXiv:0711.0798 [hep-th]]. 
[13] P. Olesen, Phys. Lett. B 265, 361 (1991) [Erratum-ibid. B 267, 541 (1991)]. 
[14] L. A. Caffarelli and Y. Yang, Commun. Math. Phys. 168, 321 (1995). 

[15] G. V. Dunne, Self-Dual Chern-Simons Theories, (Lecture Notes in Physics Monographs, Springer, Germany, 1995); G. V. 
Dunne, Aspects of Chern-Simons theory, in: Les Houches Summer School in Theoretical Physics, Session 69: Topological 
Aspects of Low-dimensional Systems, (Les Houches, France, 1998). [arxiv: hep-th/9902115]. 

[16] G. Barton and N. Dombey, Nature 311, 336 (1984). 

[17] G. Barton and N. Dombey, Ann. of Phys. 162, 231 (1985). 

[18] K. A. Milton and Y. J. Ng, Phys. Rev. D 42, 2875 (1990). 

[19] K. A. Milton and Y. J. Ng, Phys. Rev. D 46, 842 (1992). 

[20] T. -Y. Zheng, L. Zheng and C. -S. Li, Commun. Theor. Phys. 36, 431 (2001). T. -Y. Zheng, Phys. Lett. A 305, 337 (2002). 

[21] D. T. Alves, E. R. Granhen, J. F. Medeiros Neto and S. Perez, Phys. Lett. A 374, 2113 (2010). 

[22] J. Ambjorn and S. Wolfram, Annals Phys. 147, 1 (1983). 

[23] S. -S. Feng, X. -J. Qiu and Z. -Y. Zhu, Int. J. Theor. Phys. 37, 2105 (1998). 

[24] A. S. Goldhaber and M. M. Nieto, Phys. Rev. Lett. 26, 1390 (1971). N. M. Kroll, Phys. Rev. Lett. 26, 1395 (1971). 
[25] R. Jackiw and E. J. Weinberg, Phys. Rev. Lett. 64, 2234 (1990). 

[26] Sumathi Rao, "An Anyon Primer", Arxiv: hep-th/9209066; R. Banerjee and P. Mukherjee, Nucl. Phys. B478, 235 (1996); 

Prog. Theor. Phys. 101, 1189 (1999). 
[27] X. G. Wen and A. Zee, Phys. Rev. Lett. 62, 1937 (1989). 

[28] C. P. Burgess and B. P. Dolan, Phys. Rev. B 63, 155309 (2001) [hep-th/0010246]. 

[29] P. W. Anderson, Phys. Rev. 109, 1492 (1958); ibid. 130, 439 (1963). 

[30] M. Tinkham, Introduction to Superconductivity, (McGraw- Hill, New York, 1996). 

[31] A. Allocca, et ai, arXiv:1205.2478 [hep-th]. 

[32] C. M. Wilson, et. ai, Nature 479, 376 (2011). 



